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ABSTRACT 

Using eight dark matter haloes extracted from fully-self consistent cosmological 
A^-body simulations, we perform microlensing experiments. A hypothetical observer is 
placed at a distance of 8.5 kpc from the centre of the halo measuring optical depths, 
event durations and event rates towards the direction of the Large Magellanic Cloud. 
We simulate 1600 microlensing experiments for each halo. Assuming that the whole 
halo consists of MACHOs, / = 1.0, and a single MACHO mass is tom = 1-0 Mq, the 
simulations yield mean values of r = 4.7i|"2 x 10~ 7 and T = 1.6tJ'g x 10" 6 events 
star -1 yr -1 . We find that triaxiality and substructure can have major effects on the 
measured values so that r and T values of up to three times the mean ca n be found. If 
we fit our values of r and T to the MACHO collaboration observations ijAlcock et alJ 
12000(1 . we find / = 0.23±°;^ and m M = 0.44±°'^. Five out of the ei § ht haloes under 
investigation produce / and tum values mainly concentrated within these bounds. 



1 INTRODUCTION 



\ Experiments suc h as MACHO jAlcock et ail l2000tl and 
POINT-AGAPE JCalchi Novati et al.ll2005h have detected 
microlensing events towards the Large Magellanic Cloud 
(LMC) and the M31, supporting the view that some fraction 
of the Milky Way dark halo may be comprised of massive as- 
tronomical compact halo objects (MACHOs). The fraction 
of the dark halo mass in MACHOs, as well as the nature 
of these lensing objects and their individual masses can be 
constrained to some extent by combining observations of 
micr olensing events with analytical models of the dark hal o 
(e.g. lKerinslll99gl : lAlcock et al]l200d: ICardone et afll200ll) . 
However, the real dark halo might not be as simple as in an- 
alytical models, which typically have spherical or azimuthal 
symmetry and, in particular, smoothly distributed matter. 
Cosmological simulations indicate that dark matter haloes 
are neit her isotropic nor ho mogeneous but are rather triax- 
ial (e.g. IWarren et alJll992|l and conta in a notable amount 
of substructure (e.e AK\win et al.lll9 99h Even for the Milky 
Way it is still unclear if the enveloping dark matter halo has 
spherical or triaxial morphology: a detailed investigation of 
the tidal tail of the Sgr dwarf galaxy su pports the notion 
of a nearly spherical Galactic potential dlbata et alJl2001i 
iMaiewski et all2003tl whereas that the data are also claimed 
to be consistent with a prolate or oblate halo dHefmill20oj) . 
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The shape of the halo has an effect on the predicted number 
of microlenses along different lines-of-sight, for comparison 
with the results of microlensing experiments. 

In this paper, we examine the effects of dark halo mor- 
phology and dark halo dumpiness on microlensing surveys 
conducted by hypothetical observers located in eight N- 
body dark matter haloes formed fully self-consistently in 
cosmological simulations of the concordance ACDM model. 
The eight haloes are scaled to approximately match the 
Milky Way's dark halo in mass and rotation curve prop- 
erties, and hypothetical observers are placed on the surface 
of a "Solar sphere" (i.e. 8.5 kpc from the dark halo's centre) 
from where they conduct microlensing experiments along 
lines-of-sight simulating the Sun's line-of-sight to the LMC. 
The effects on microlensing of triaxiality and clumping in 
the haloes are examined. 

The paper expands upon work in a study by 
IWidrow fc Dubinskil lll998f) although we take a slightly dif- 
ferent point of view. Instead of studying the errors caused 
by different analytical models fitted to an Af-body halo, we 
analyse the microlensing survey properties of the simulated 
haloes directly. That is, we do not try to gauge the credibil- 
ity of analytical microlensing descriptions but rather use our 
self-consistent haloes for the inverse problem in microlensing 
JCardone et all 1200 ill . In other words, the free parameters 
(the MACHO mass and the fraction of matter in MACHOs) 
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are determined from the observations via the optical depth, 
event duration and event rate predictions of the models. This 
not only allows us to make predictions for these parameters 
based upon fully self-consistent halo models but also to in- 
vestigate the importance of shape and substructure content 
of the haloes. 

The paper is structured as follows. We describe the N- 
body haloes and their preparation in Section 2, the sim- 
ulation details are covered in Section 3, equations for the 
observables are derived in Section 4, results are given in 
Section 5, some discussion is presented in Section 6, final 
conclusions are in Section 7. 



2 THE HALOES 

2.1 Cosmological simulation details 

Our analysis is base d on a suite of eight high-re solution 
iV-body simulations jGill. Knebe fc Gibson! l2004al) carried 
out using th e publicly available adaptive m esh refinement 
code MLAPM (Kncbc. Green fc BinnevlEoOll) in a standard 
ACDM cosmology (fi = 0.3, fl x = 0.7, Q b h 2 = 0.04, h = 
0.7,(78 = 0.9). Each run focuses on the formation and evo- 
lution of galaxy cluster sized object containing of order one 
million collisionless, dark matter particles, with mass reso- 
lution 1.6 x 10 s hT 1 Mq and force resolution ~2/i -1 kpc (of 
order 0.05% of the host's virial radius). The simulations have 
sufficient resolution to follow the orbits of satellites within 
the very central regions of the host potential 5-10 % of 
the virial radius) and the time resolution to resolve the satel- 
lite orbits with good accuracy (snapshots are stored with a 
temporal spacing of At ~ 170 Myr). Such temporal resolu- 
tion provides of order 10-20 timesteps per orbit per satellite 
galaxy, thus allowing these simulations to be used in a pre- 
vious paper to accurately measu re the orbital par ameters of 
each individual satellite galaxy jGill et al.ll2004bT) . 

The clusters were chosen to sample a variety of en- 
vironments. We define the virial radius i? v ir as the point 
where the density of the host (measured in terms of the 
cosmological background density pt) drops below the virial 
overdensity A v i r = 340. This choice for A v i r is based upon 
the dissipationless spherical top-hat collapse model and is 
a function of both cosmological model and time. We fur- 
ther applied a lower mass cut for all the satellite galax- 
ies of 1.6 x 10 10 h' 1 M Q (100 particles). Further specific 
details of the host haloes, such as masses, density pro- 
files, triaxial ities, environment and merge r hist ories, can 
be foun d in iGill. Knebe fe Gibson! J2004a!) and iGill et alJ 
J2004bl) . Table gives a summary of a number of rele- 
vant global properties of our halo sample. There is a promi- 
nent spread not only in the number of satellites with mass 
M sat > 1.6 x 10 10 /i _1 M Q but also in age, reflecting the dif- 
ferent dynamical states of the systems under consideration. 



2.2 Downscaling 

As we intend to perform microlensing experiments directly 
comparable to the results of the MACHO collaboration we 
need to scale our halo es to the size of the Milky Way . For 
this purpose we follow iHelmi. White fc Springel (120031) and 



Table 1. Summary of the eight host dark matter haloes. The 
superscript cl indicates that the values are for the unsealed clus- 
ters. Column 2 shows the virial radius, R v i T ; Column 3 the virial 
mass, M v ; r ; Column 4 the redshift of formation, Zf orm ; Column 5 
the age in Gyr; and the final column an estimate of the number 
of satellites (subhaloes) in each halo, iV sa t. 



Halo 


E? cl 
Jtvir 


M ■ cl 

iVI vlr 


form 


age 


AT sat 




h' 1 Mpc 


10 14 h- 1 Mq 




Gyr 




# 1 


1.34 


2.87 


1.16 


8.30 


158 


#2 


1.06 


1.42 


0.96 


7.55 


63 


#3 


1.08 


1.48 


0.87 


7.16 


87 


#4 


0.98 


1.10 


0.85 


7.07 


57 


# 5 


1.35 


2.91 


0.65 


6.01 


175 


#6 


1.05 


1.37 


0.65 


6.01 


85 


# 7 


1.01 


1.21 


0.43 


4.52 


59 


#8 


1.38 


3.08 


0.30 


3.42 


251 



apply an adjustment to the length scale, requiring that den- 
sities remain unchanged. Hence the scaling relations are: 
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where r is any distance, v is velocity, m is mass, t is time 
and the superscript cl refers to the unsealed value. Because 
the haloes are downscaled, 7 will always be in the range 
7 6 [0,1]. 

To find a suitable 7, the maximum circular velocity of 
each halo is required to be 220 kms -1 . The resulting (scaled) 
rotation curves for all eight haloes are shown in Figure 
This figure highlights that Halo #8 is a special system — its 
evoluti on is dominated b y the interaction of three merging 
haloes jGill et alJl2004bl) . We do not consider Halo #8 to 
be an acceptable model of the Milky Way, but we include it 
into the analysis as an extreme case. 

One might argue that our (initially) cluster sized ob- 
jects should not be used as models for the Milky Way for 
they formed in a different kind of environment with less 
time to settle to (dynamical) equilibrium (i.e. the oldest of 
our syste ms is 8.3 Gyr vs. ~ 12 Gyr for the Milky Way). 
However. lHelmi. White fc Springel j200St) showed that more 
than 90% of the total mass in the central region of a realis- 
tic Milky Way model was in place about 1.5 Gyr after the 
formation of the object. Our study focuses exclusively on 
this central region (i.e. the inner 15% in radius) boosting 
our confidence that our haloes do serve as credible models 
of the Milky Way for our purposes. 

In Table|2]we list the scaled radii and masses, as well as 
the scaling factor 7, for each halo. Note that we have adopted 
a Hubble constant of h = 0.7, which applies throughout 
the paper. The circular velocity (rotation) curves for the 
scaled haloes are shown in Figure and demonstrate how 
the scaled mass M(< r) is accumulated out to 400 kpc. Two 
vertical lines mark the position of observers at the Solar 
circle and the distance of the LMC from the halo centre; 
within this region, we note that the mass profiles are similar 
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Table 2. Properties of the haloes after downscaling to Milky Way 
size. Column 2 shows the scale factor for each halo, 7; Column 
3 shows the scaled virial radius, i? v i r ; Column 4 the scaled virial 
mass, M v j r ; and the final column the scaled particle mass, m v . 



Halo 


7 




M vir 


m p 






kpc 


10 12 M 


10 6 M 


#1 


0.197 


380 


3.08 


1.75 


#2 


0.248 


378 


3.00 


3.47 


#3 


0.253 


391 


3.39 


3.71 


#■1 


0.277 


388 


3.23 


4.86 


#5 


0.197 


383 


3.14 


1.76 


#6 


0.267 


402 


3.61 


4.33 


#7 


0.278 


403 


3.67 


4.92 


#8 


0.213 


421 


4.19 


2.21 



240 




r (kpc) 

Figure 1. Circular velocity curves of the downscaled haloes 
(V c = \J GM(r)/r). Note that the maximum of each curve is 
at 220 kms -1 . The scaling factors 7 are shown in Table EJ The 
two vertical lines at 8.5 kpc and 50.1 kpc mark the distances of 
the Sun and the LMC from the centre of the Milky Way. 

in all haloes, however, there are differences, especially for 
halo #8 (the dynamically merging system of haloes) . 

2.3 Characteristics 

Dark matter haloes within cosmological simulations are very 
different from the idealised analytical dark haloes that are 
described by only a few parameters. Two significant char- 
acteristics of the dark haloes that are not captured by the 
analytical models are triaxialty and substructure. For ex- 
ample, within the simulations on both cluster and galac- 
tic scales, the dark matter haloes are triaxial r ather than 
spherical dWarren et alJll992t Ijing fc Sutoll2002f) . However, 
when gas dynamics is included within the simulation, these 
clusters /galaxies a r e more spherica l IIKatz fc GunrJ Il99ll 
lEvrard et all Il994 iDubinskil IT99I . and when cooling is 
included within the simulation, t he haloes are even more 
spherical jKazantzidis et aT1l2004l) . 

Simulations show that CDM haloes contain a large 
number of sub haloes at both ga l actic and cluster scales 
(i.e. satellites) jKlvpin et alJll999l fMoore et alJll999lh fur- 
ther lMoore et alTiuggsF " assert that galaxy haloes appear as 



scaled versions of galaxy clusters with respects to the sub- 
halo populations. 

However, within the inner regions of the host haloes (on 
both galactic and cluster scal es) the number density of sub- 
haloes decreases considerably <|G higna et al.l2000LlGill et alJ 
l2004bl|piemand. Moore, fc Stadel2004) . Even though these 
results seem to hold under num erical convergence tests 
JPiemand. Moore, fc Stadel2004) . the dynamical properties 
of the dark matter subhaloes differ from the observed galax- 
ies within clusters. 

With the use of semi-analytical techniques these 
differences have bee n accounted for iGao et al.l 120041 : 
iTavlor fc Babulll2004l) . iGao et al.l d2004l) conclude that the 
tidal stripping of the dark haloes was very efficient in reduc- 
ing the dark matter mass of the subhaloes while having less 
of an effect on the tightly bound baryonic galaxy at its core. 
Thus, while the dark matter subhaloes in the inner regions 
were being disrupted, baryonic cores surv ived. 

For our pure dark matter simulations jGill et all2004bft 
we also see this lack of dark matter substructure in the inner 
regions. This is shown in Table|3] where we list the subhaloes 
(along with some of their (downscaled) integral properties) 
found in the inner 70 kpc of each halo. Note that masses 
of the subhaloes are limited to a relevant range, 5.0 x 10 s 
M < M sat < TO x 10 10 M . 

For the Milky Way, we know that at least one large 
satellite exists within the inner 70 kpc: the LMC. Thus, 
whatever the paucity of satellites in inner regions of 
dark haloes in the simulations, here at least is one real- 
world example. The mass of the LMC is estimated to be 
M LM c(8.9 kpc) = 8.7T4.3 x 10 9 M from kinematical data 
jMarel et alJl2002l) . Reassuringly enough, we do find a few 
satellites with similar masses (cf. Table |3J in the simulated 
haloes. Whether there are more (dark) satellites in the in- 
ner regions of the Milky Way than in our models is an open 
question. 

3 SETTING THE STAGE 
3.1 The microlensing cone 

A Galactocentric coordinate system is given by the three 
eigenvectors of the inertia tensor of the respective halo. We 
choose the z-axis to coincide with the major axis. To con- 
duct our microlensing experiments, we then place an ob- 
server at the Solar distance from the centre of a halo. There 
are various options for possible orientations with respect to 
the coordinate system; we basically restrict ourselves to the 
scenario where the observer is placed somewhere on a sphere 
of radius 8.5 kpc. 

We sample the halo with a number of sightlines and 
combine the individual results obtained from these indepen- 
dent measurements to obtain values of interest. Once we put 
down an observer, we can construct an observational cone 
towards a region simulating observations of an LMC sized 
patch of sky, and compute the microlensing signal generated 
by particles in the cone out to a distance of 50.1 kpc (see 
Figure |5J. The hypothetical LMC has also to be located in 
a realistic direction in respect to the observer. 

We use 40 observers uniformly distributed on a sphere 
of radius 8.5 kpc. Each observer observes 40 uniformly sepa- 
rated LMCs, for which the angle (r b s , tlmc — r h s ) is kept 
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Table 3. Substructure within the inner 70 kpc of the (down- 
scaled) haloes. Column 1 shows the name of the halo; Column 2 
the distance from the centre of the subhalo to the centre of the 
host, Z sa t; Column 3 the truncation radius of the subhalo (within 
the truncation radius, the density of the subhalo is larger than the 
background density of the host halo) , r sa t ; Column 4 the number 
of bound particles in the subhalo, JV sa t; and the final column the 
mass of the subhalo, M s &t- 



Halo 


Zsat 


^sat 


iVsat 






kpc 


kpc 


particles 


1O 9 M 


a -i 
#J- 


/ .U 


Q/l P. 

O4.0 




o on 


#1 


20.0 


45.4 


2957 


5.27 


#1 


23.5 


27.9 


632 


1.13 


#2 


53.5 


24.1 


207 


0.73 


#2 


44.8 


23.3 


186 


0.66 


#3 


60.7 


38.1 


793 


2.99 


#3 


63.0 


32.7 


466 


1.76 


#4 


67.8 


30.8 


317 


1.57 


#7 


38.5 


48.7 


1217 


6.09 


#7 


62.4 


22.6 


124 


0.62 


#8 


46.1 


32.9 


859 


1.93 



fixed at 89.1° (the original angle between (f©, ?l,mc—Tq) in 
Galactocentric coordinates). This configuration (illustrated 
in Figure gives 1600 individual sightlines, and microlens- 
ing cones, which sample different sets of particles within a 
halo. The choice for the number of sightlines is a compro- 
mise between computation time and sampling density. As 
seen in Figure [3] 1600 sightlines sample the volume quite 
sufficiently, even when the sightlines are drawn as lines in- 
stead of actual cones. 

The microlensing cone itself is defined by the observer 
and the source. Figure [5] illustrates that we treat the source 
(i.e. LMC) as a circular disk with diameter ds- The parti- 
cles inside the cone are considered to be gravitational lenses. 
However, this kind of setup can lead to a handful of parti- 
cles (those closest to the observer) contributing most to the 
optical depth. This leads to high sampling noise, and is es- 
sentially due to the limited resolution of the cosmological 
simulation (see Figure 2J. It follows that the particle distri- 
bution has to be smoothed in some manner to suppress this 
noise; we describe this and the related issue of the MACHO 
masses in the next section. 



3.2 Mass resolution vs. MACHO mass 

In order to compute microlensing optical depths and other 
properties from the particles in the simulations, we need 
to get from a mass scale of order 10 6 Mq, i.e. that of the 
particles in the simulation, down to the mass scale of the 
MACHOS, of order 1 Mq; i.e. we account for the fact that 
a typical particle in the simulations represents of order 10 6 
MACHOs. 

To overcome the limitation imposed by the mass reso- 
lution of the simulation we recall that individual particles in 



Source with an solid angle of 

(o s = n(r s /lf 



The microlensing cone 




Observer 



Particle with an solid Einsten 
angle of co = n(r E /l) 2 



Figure 2. The microlensing sightline is defined by the locations 
of the observer and the source and the cone by the sightline and 
the source diameter. We use Zg = 50.1 kpc and dg = 4.52 kpc. 
The particles inside the cone are treated as gravitational lenses. 
Typically ~ 100 particles from the cosmological simulation are 
found inside a cone. This number is increased to ~ 10000 by 
breaking them up into subparticles according to the recipe out- 
lined in Section l3.2l 



the simulation are not treated as 5-functions but have a fi- 
nite size. The extent of a particle (i.e. its size) is determined 
— in our case — by the spacing of the grid, as we are us- 
ing an adaptive mesh refinement code (i.e. MLAPM). In MLAPM 
the mass of each particle is assigned to the grid via the 
so-calle d triangular-shaped cloud (TSC) mass- assignment 
scheme (iHocknev fc Eastwoodll98lT) which spreads every in- 
dividual particle mass over the host and surrounding 3 3 — 1 
cells. The corresponding particle shape (in ID) reads as fol- 
lows: 



S(x) = \ 1 ~ ~ im \ x \ <L 

otherwise 



(5) 



where L is the spacing of the grid and x measures the 
distance to the centre of the cell. 

For each particle present in and around a particular 
cone along an observer's line-of-sight, we determine the size 
of the finest grid surrounding it, which in turn determines 
the physical extent of the particle. The particle is then re- 
sampled with 100 "subparticles" whose positions are ran- 
domly chosen under the density distribution given by Equa- 
tion |S] Figure |5] illustrates these "subparticle clouds", sam- 
pled according to this TSC mass- assignment. 

The improvement gained by using this method can be 
viewed in Figure 0] where we show measured optical depths 
(to be defined in Section 14.21 in a series of microlensing 
cones towards the LMC for observers continuously rotated 
on the edge of a disc perpendicular to the major axis of the 
dark halo. What is shown here is the noisy "original" esti- 
mates of optical depth for many lines-of-sight, compared to 
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Figure 3. Sightlines used in the microlensing simulation. The 
image shows all 1600 sightlines of the 40 observers. The observers 
are located on a sphere with a radius of r© = 8.5 kpc. The sphere 
can be seen as the "dark ball" in the centre of the sightlines. The 
radius of the sphere which holds the sources is 49.5 kpc. Every 
observer has 40 sources (LMCs) at a distance of 50.1 kpc. The 
sightlines are not uniformly distributed because there is a limited 
number of observers and the angle (r b s , r^MC ~ fobs) ls kept 
fixed at 89.1°. This is the original angle between (tq, rLMC - ^) 
in Galactocentric coordinates. 



the optical depth measurements obtained when the particles 
within the cones have been resampled to "subparticles" . In 
the un-resampled case, the handful of particles which hap- 
pen to reside close to the observer are found to dominate 
the calculation of the optical depth and other quantities of 
interest such as event duration and event rate. Variations in 
the optical depth from sightline to sightline can vary by up 
to a factor of three in the un-resampled sample cones, sim- 
ply because a handful of particles dominate the microlens- 
ing. Convergence experiments have shown that breaking the 
original simulation particles down to 100 subparticles is suf- 
ficient to reduce the noise from this source to be negligible; 
if we were to use more subparticles, there is no further gain 
in accuracy. 

The subparticles introduced have masses of order 10 4 
M© , so we are unfortunately still orders of magnitudes away 
from actual MACHO masses. In the formulae introduced 
later (Section QJl, each subparticle is expressed as a sin- 
gular concentration of MACHOs. That is, a subparticle is 
treated as a set of MACHOs which have the same posi- 
tion and velocity as the subparticle itself, and for conve- 
nience in the simulations the MACHO mass is chosen to be 
Tom = 1.0 Mq (although this is later relaxed by scaling). Fig- 
ure demonstrates the hierarchy of particles starting from 
the cosmological simulation particle down to the individual 
MACHOs. 
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Figure 4. Comparison between the use of cosmological simula- 
tion particles and triangular-shaped clouds containing 100 sub- 
particles. For the illustration, we distributed 132 observers on a 
perimeter of a disk with a radius of 8.5 kpc, located on the xy- 
plane. The volumes of consecutive cones overlap by half. Without 
any resolution enhancements the variation in optical depth from 
cone to cone can be as large as a factor of three! This is due 
to the insufficient mass resolution of the cosmological simulation. 
We reached the Poissonian noise level, in which cone-to-cone vari- 
ations are < 20 %, with 100 subparticles, by experimenting. The 
large scale trend is due to triaxiality of the halo. 



4 THE MICROLENSING EXPERIMENTS 
4.1 Definitions 

Before going to the microlensing equations, we define what 
we mean by certain terms: 

A microlensing source is a circular region, oriented 
perpendicularly to the line-of-sight of the observer, in which 
uniformly distributed background source stars are located. 
Because of the statistical nature of the source model, the 
number of background stars affects only the microlensing 
event rate. The region with a diameter ds and a distance 
Is, is shown in Figure [5] 

A source star is a star located somewhere in the disk 
of the microlensing source. 

A microlens is also a circular region, oriented perpen- 
dicularly to the line-of-sight of the observer and centered 
on a dark particle. The region has a radius te, which de- 
pends on the location and mass of the particle. Lenses are 
always inside a microlensing cone between the source and 
the observer. 

A microlensing event is a detectable amplification of 
a source star caused by a microlens. Detectable means that 
the light of a source star is amplified by a factor larger than 
1.34. This occurs when the sightline to a source star passes 
a particle within the lens radius, also known as the Einstein 
radius. 



4.2 Optical depth 

The Einstein radius rE of a gravitational lens is defined as 

■*-^(i-s)-^(i-E> 
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1 .2 kpc box 




Figure 5. For this illustrative example of the subparticle clouds, 
we placed 222 cosmological simulation particles onto the xj/-planc 
and broke each particle down to 10,000 subparticles. To get the 
subparticle positions, we used the triangular-shaped cloud den- 
sity profile identically to the actual cosmological simulation. The 
subparticles form a cube (of a size of the corresponding MLAPM 
grid cell) around the position of the cosmological simulation par- 
ticle. In the real microlensing simulation, the whole volume under 
investigation is filled with intersecting cubes and the void areas 
between the cubes seen here are not present. 
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PARTICLES 
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BREAKING UP TO MACHOs 
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m, 
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Figure 6. Splitting the cosmological simulation particle to 100 
subparticles and a subparticle to JVm MACHOs. In the first phase, 
we break the cosmological simulation particle to 100 subparticles. 
The resulting subparticle represents an order of JVm ~ 10 4 MA- 
CHOs with a mass of 1 Mq . The subparticle is treated as a MA- 
CHO concentration where the represented MACHOs all have a 
single location (£m = ^sub)i mass (myi = m^-^/Nyi) and velocity 
(»M = "sub)- 



where m is the mass of the lens, I is the distance to the 
lens, Is is the distance to the source and A = 

For a single subparticle, we define the lens radius as 



^sub'sub 



^sub 



(7) 



where m su b is the mass of the subparticle and Z su b the 
distance between the subpartice and the observer. Note that 
all subparticles within a halo have the same mass. 

The optical depth of a single lens describes the proba- 
bility that any given source star is amplified by the lens at 
any given time. Thus, the general equation for the optical 
depth of a lens can be written as the ratio of the solid angles 
of the lens and the source 



(8) 



where u> is the solid Einstein angle of the lens and u>s 
is the solid angle of the source. When the solid angles are 
small, this can be approximated to 



rg Is 
I r s 



AmO~ 



(9) 



where rs is the radius of the source and 0~ 2 = (ls/ r s) 2 - 
Note that is approximately the half opening angle of the 
cone. 

The optical depth of a subparticle with a mass m su b is 
defined as 



",1,1. - 10 2 m sub [ 

'•sub "8 



(10) 



We take this value to be the total optical depth from 
all the ~ 10 4 MACHOs the subparticle represents. 

Optical depth is additive, as long as the lenses do not 
overlap or cover the whole source, and so the total optical 
depth in a microlensing cone which contains N subparticles 



N 



'sub 



= AO 



N 



(11) 



From Eauation llll it follows that the optical depth in a 
cone depends mainly on the distances of the subparticles re- 
lated to the observer. The closer the particles are the larger 
is the optical depth. Section 13.21 covered the details of par- 
ticle breaking, which also has a large effect on r CO no. Note 
that T conc does not depend on the MACHO mass Tom- 

4.3 Event duration 

The event duration of a lens describes the typical duration of 
the amplifying event the lens would produce. The detected 
event durations in the MACHO experiment are the order of 
100 days. Event duration depends on the tangential velocity 
with which a lens would seem to pass a source star and on the 
Einstein radius of the lens. The equation for an individual 
subparticle is 

tsub = TT - ^! (I 2 ) 
2 v Buh 

where r su b is the lens radius of the subparticle and v an h 
the apparent tangential velocity difference between the sub- 
particle and the source, respective to the observer. 
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The term ^ is the average crossing length of a circle 
with an unit radius. By adding this term, we correct the 
event duration from the maximum value (2rE/v) to the sta- 
tistically expected value. The correction can also be seen in 
the observational optical depth, r mea3 , as a 7r/4 term in Eq. 
1 of Alcock et al. (2000). 

Using Equation |7| Equation 1121 becomes 



tsub 



1/2 ^sub 
Usuh 



.4 



1 



1 1/2 



(13) 



However, this is not yet the value we are after. We need 
to solve tsub, which is the average event duration caused by 
the subparticle's MACHOs. To get t S ub, we simply replace 
m su b with mm and assign Im = £sub and «m = «sub- We get 



, _ 7T 1/2 ^sub 

t M = ~m M - — 

2 Vsub 



A 



1 

^sub 



1/2 



(14) 



where tM is the event duration of a single MACHO when 
the MACHO inherits the location and the velocity of the 
subparticle. What is essentially stated in Equation ll4l is that 
we use the event duration of a single MACHO as the average 
event duration for the whole subparticle. 

As is well known, Equation 1131 shows that, unlike the 
optical depth, event duration is a function of the MACHO 
mass. This fact can be used to find a preferred MACHO 
mass for a given event duration. 



4.4 Event rate 

Event rate F for a given experiment is simply the number of 
expected events per an observing period for a given source. 
For a single lens it can be given as T — r/t, i.e., as the ratio 
of the optical depth and event duration. 

The MACHO collaboration observed ~ 13 - 17 events 
in 5.7 years. Usually, T is given in events star -1 yr _1 . The 
MACHO collaboration had 10.7 x 10 6 source stars which 
gives (a maximum) of F = 17/5.7/(10.7 x 10 6 ) events star -1 
yr -1 = 2.79 x 10 -7 events star -1 yr -1 . 

We calculate the event rate for a subparticle consisting 
of N M MACHOs as 



r AT T AT ™ AT T sub/N M T sub 

1 sub = "Ml M = JVM- — = J»m- 

IM 



^sub 



tsub 



(15) 



where Tm = p* 1 is the event rate for a single MACHO. 
By using Eouations' ll0l and ll3l the last term of Eouation ll5l 

expands to 



2 

7T 



i ~ /)1/2q-2 -1/2 "sub 

Tsub = -A O m M m su b 



1 1 

^sub \ ^sub Is 



1/2 



For the cone, we simply sum the event rates of the sub- 
particles together, as was done to the optical depths in Equa- 
tion [TTl 



-A 1/2 e" 2 r 



-1/2 



TTlsub 



N (i) , \ 1/2 

V^^sub/ 1 1 



sub x sub 



Note that r su b oc m^, 1 ' 2 whereas t S ub ex m^ 2 . 



5 RESULTS 

5.1 Differential optical depths and event rates 

The results of our simulations are compared to the obser- 
vations via the event duration because this is the only di- 
rectly observed quantity. Our aim is to determine how optical 
depth r and event rate F depend on event duration. 

In Figure[5] we show dr/dt and dT/dt, averaged over all 
the 1600 cones for each halo, i.e we show for each simulated 
halo the global behaviour of differential optical depths and 
event rates as a function of event duration. The normalised 
standard deviations in these quantities amongst sightlines 
are shown (in separate panels, for clarity). Both curves for 
r and F peak in the range 50 and 100 days. This was ex- 
pect ed as it is t he result obtained from analytical models 
(e.g. lKerinsll998() . Our simulations and the analytical results 
are in this sense adequate fits to the actual observations, in 
which the (few) events cluster around 100 days duration. 

In Tablc0] we list the modal values for the distributions, 
t cxp , and the averaged cone values, (r conc ) and (r c0 nc), which 
are numerically integrated values of the differential functions 
in Figure |5] Note that the error limits on the last row (la- 
belled "mean") are simple mean values of the error limits 
from the corresponding column because they represent the 
mean variations between cones in a halo, not variations be- 
tween haloes. 

It can be seen that the integrated values differ between 
haloes. One might argue that there is a trend towards lower 
values from Halo #1 to Halo #8, and this is perhaps no 
surprise because the haloes were ordered by their dynamical 
age. Similar grouping of curves can be seen in the circular ve- 
locity profiles in Figure Q Because circular velocity profiles 
reflect density profiles, it is the differences in matter densi- 
ties between haloes that can be seen to cause differences in 
the integrated values of dr/dt and dT/dt. 

We have applied the MACHO collaboration efficiency 
functions to our simulations in order to get an idea of how 
their efficiency affects the overall observables. The results 
can be seen in Tables|S]and|n] The efficiency functions reduce 
(Tconc) and (r cono ) by more than a factor of two and shift 
the expected event durations to be longer. 

Note that all the following results are given in terms 
MACHO's efficiency function "A" — this is the more con- 
servative option of two efficiency functions given by the MA- 
CHO collaboration; tests show that both give very similar 
results. 



5.2 Substructure 

An interesting feature can be seen in the relative standard 
(lg) deviation (i.e. normalised to the global mean) in Figure[5]of 
both dr/dt and dT/dt for Halo #1. The relative standard 
deviation of both functions has a peak around t — 50 days of 
greater than 1.0 compared to typical values of ~ 0.7. This 
peak is due to a subhalo with a mass of 5.27 x 10 9 M , 
i.e. the second subhalo of Halo #1 listed in Tabled] It is 
the only satellite with a substantial mass, concentration and 
(17) location to cause any features to the r and T distributions. 

For example, Halo #7 contains a subhalo with a comparable 
mass (6.09 x 10 9 M ), but because it is located further away 
from the centre of the halo (and from the observer), it fails 
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Table 4. The mean optical depths and event rates, with their 
standard deviations, over the 1600 cones for each halo. See Figure 
[TJfor the corresponding (differential) distribution functions — the 
numbers represented here are the integrated values of those func- 
tions. We also list the modes of the distributions: i cxp columns 
give the most probable event duration for the respective distri- 
bution. Error limits show the variations between cones (i.e. mi- 
crolensing experiments) in a halo. 



Halo 


(Tcone ) 


^cxp 


(Tccmc } 


^exp 




io- 7 


days 


10~ 6 


days 








events 










star -1 yr — 1 





#1 

#2 
#3 
#■1 
#-> 
#6 
#7 
#8 



5.6 
5.9 
1.3 
5.1 
4.1 
5.0 
4.7 
2 



6.0 
2.7 
+7.5 
2.7 
+3.8 
1.8 
+5.6 
2.4 
+3.4 
2.2 
+4.9 
1.9 
(i.fi 
2.5 
+2.1 
-1.0 



84 
74 
84 
86 
91 
81 
84 
104 



2.0 
2.1 
1.4 
1.6 
1.3 
1.7 
1.6 
0.8 



+1-7 
0.8 
+2.1 
0.8 
+1.1 
0.5 
+ 1.3 
0.7 
+0.9 
0.6 
+ 1.3 
0.6 
+ 1.6 
0.7 
+0.6 
0.3 



61 
61 
64 
69 
71 
61 
81 



4.7 



+5.0 



86 



1.6 



+ 1.3 



to affect the overall r and T distributions. The subhalo in 
Halo #1 is therefore the only substructure which could be 
directly associated with microlensing signal not associated 
with the global properties of the haloes. 

It should be emphasized that the distributions shown 
in Figure [5] are averaged over all cones. The subhaloes, if 
present in a halo, contribute only to a small number of the 
cones and thus, possibly peculiar distributions are lost by 
the averaging process. In another words, the distribution 
functions of a single cone can be quite different from the 
averaged one. A particularly severe case of this can be seen 
in Figure[7]where we show the average distribution functions 
of 24 cones which penetrate the second subhalo in Halo #1 
compared to the overall distributions. The cones are chosen 
so that their lines-of-sight are closer than 5 kpc to the centre 
of the subhalo and therefore the cones contain particles from 
the subhalo 's core. The cones yield r and F values which are 
approximately twice as large as the average values in Table 
0] Moreover, the event duration peaks (t exp ) are different 
from the Table[I]values because the subhalo's particles have 
different velocities than the particles of the host halo. 

We also probed the cones which penetrate the first sub- 
halo in Halo #7 for signatures in the differential functions. 
As mentioned earlier, this subhalo has a comparable mass 
to the major subhalo of Halo #1. We found that an excess 
signal is also present in the "subhalo cones" of Halo #7. 
However, this subhalo does not produce such a peaked func- 
tions as seen in Figure[7|but instead, dr/dt and dT/dt stay 
steadily above the average shape, at two times the average 
value, between 40 days < t < 200 days. This roughly dou- 
bles the values of r and F for cones penetrating the subhalo, 
compared to the average cone values. 




100 120 140 160 180 200 
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20 
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Figure 7. "Subhalo cones" shows the average differential func- 
tions from 24 cones which penetrate the second subhalo in Halo 
#1 (see Tablel3"l. "All cones" is shown for comparison — it is the 
the same function as in Figure|5]for Halo #1. The large deviation 
from the average shape is seen in the standard deviation panels 
of Figure [5] and here as a sharp rise of both dr / dt and dr / dt 
functions starting from t ~ 40 days. Individual "subhalo cones" 
produce even more distorted distributions and we had to average 
over the 24 cones to get the peak to stand out clearly. 



5.3 Triaxiality 

In addition to substructure, triaxiality was the other halo 
component which we wanted to study from the microlens- 
ing point of view. Triaxiality effects are not seen in the dr/dt 
and dT/dt functions because they are averaged over all the 
cones, and the cones are distributed spherically inside a halo. 
However, when the cones are not distributed spherically, one 
can immediately see effects due to triaxial shape. Figure 2] 
shows how an observer on a "Solar orbit" observes a roughly 
sinusoidal variation for the optical depth as function of posi- 
tion on the perimeter. We find the amplitude to range from 
T"mm ~ 4 x 10~ 7 to T max ~ 6 x 10 -7 and a peak-to-peak 
wavelength of ~ 60 cones (which corresponds to 180 degrees 
in observer positions) . This shape is solely due to triaxiality 
and is seen in all the haloes. 

To measure triaxiality related variations more quanti- 
tatively we constructed Figures 1101 and 1111 In both these 
figures, we have measured the angle between a positive co- 
ordinate axis and the source, keeping in mind that the co- 
ordinate axes are aligned with the triaxial axes. After this 
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Figure 8. An example which shows the anticorrelation between 
apparent tangential velocities and locations of subparticlcs inside 
a microlensing cone. This anticorrelation explains why triaxiality 
does not show up in Y as clearly as in r. The size range of the 
triangular-shaped clouds also shows up clearly in the plot. 



is done for all cones, we bin the observables and calculate 
the upper and lower fractiles so that they hold 65 % of the 
values. Especially the optical depth fractiles reveal clear tri- 
axiality signal. 

In the first panel in Figure [TU1 the mean optical depth 
is largest when the sources are close to the z-axis (i.e. the 
major triaxial axis), as one would expect. For haloes #1, #2 
and #7, the optical depth on opposite sides of the :ry-plane 
differs somewhat, demonstrating that matter is not neces- 
sarily distributed with azimuthal symmetry in the simulated 
(or real) haloes. 

There are high mean optical depth values near the pos- 
itive y-axis for Halo #1 in the second panel. This anomaly 
is due to the subhalo discussed in previous Section f5. 21 Now 
we can see that the subhalo is located close to the positive 
j/-axis. This is confirmed by the true location of the subhalo, 
which is found to be Z sa t = (—3.37, 19.8, 1.36) kpc. 

The lowest panel is different from the two previous ones 
because we are measuring triaxiality. The optical depth val- 
ues from sources near the x-axis seem to be the lowest 
amongst all panels. These low values (r cone ~ 4 x 10 -7 ) 
can also be seen in the second panel as the lower fractiles at 
a ~ 90°. Thus, by these "observations" we can verify that 
the z-axis is the minor axis. In the lowest panel, the fractiles 
bend upwards at a ~ 90° because the sources near both the 
z- and y-axis produce larger values than the sources near 
the i-axis. 

Triaxiality can not be seen as clearly in Figure [TT1 as in 
Figure [T01 because F is affected by the lens velocities whereas 
r is not. Apparently, the velocities do not carry enough infor- 
mation about the triaxiality and the correlation between the 
source location and the observables are somewhat "washed 
out" . This is confirmed in Figure|H|where we show that there 
is virtually no correlation between the apparent tangential 
velocity and the location of the subparticles inside a cone. 
Nevertheless, the most striking features, i.e. the 2-axis ma- 
jority and the subhalo in Halo #1, are still clearly visible 
in r. For example, a suitable subhalo can produce roughly 
twice as many events as a similar area without one, based 
on the second panel in Figure ITT1 



5.4 Estimating MACHO mass and halo mass 
fraction in MACHOs 

The MACHO collaboration used several analytical halo 
models to estimate the typical, individual MACHO mass, 
rriM, and the mass fraction of MACHOs in the halo, /, re- 
sponsible for the lensing events observed. For example, one 
of their models fits m M = 0.6±g;^ M and / = 0.21±°;J°. 
We calculate our own estimates for tum and / to see how 
our Af-body halo models compare to the analytical ones in 
making these predictions. 

As the observational constraint, we use the blending 
corrected event durations chosen by the MACHO collab- 
oration's criterion "A", t st (A), from Table 8 in Alcock et 
al. (2000). These values are not corrected for the observa- 
tional efficiency function £ (!) nor have they been reduced 
for the events caused by the known stellar foreground pop- 
ulations, given in Table 12 in Alcock et al. (2000). We 
adopt N e xp(A) = 2.67 for known population events and ar- 
rive at the observational values of r b s = 3.38 x 10~ 8 and 
Tobs = 1-69 x 10~ 7 events star -1 yr _1 . We then compare 
these values with our simulation produced values, r cone and 
Tconc, to get myi and /. The simulation values are corrected 
with the efficiency function £(t) prior to calculating the pre- 
dicted values. The calculations are based on the following 
equations. 

First, we assume that 



/ oc T, 



(18) 



which simply describes that that the probability of ob- 
serving an event is proportional to the number of lenses in 
the whole halo. From Equation 1171 and from the fact that 
the number of observed events also follows the total number 
of lenses, we get 1 



T oc fm- 



-1/2 



(19) 



These two equations permit us to estimate myi and /. 
From Equation 1181 we get 



fob 



7~obs 



(20) 



For the microlensing simulations we assumed / CO ne 
1.0 and thus, 



f=fo 



7"obs 



(21) 



This estimate of / is calculated for each cone. 
For mM, we need to use the event rate in addition to 
the optical depth. From Equations 1181 and 1191 we get 



, cone 
'M 



fob 



cone r bs 



(22) 



where / conc = 1.0 and m^ ne = 1.0 M . Thus, we get 
the equation for the MACHO mass, 



obs 



fot 



Tconc 



M . 



(23) 



1 Note that in Equation 1171 we assume that the whole halo con- 
sists of MACHOs, i.e., / = 1.0. 
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Table 5. Same as Table |1] but here we used the MACHO col- 
laboration's efficiency function A. These numbers should corre- 
spond to the E (t)— uncorrected observations if the MACHO frac- 
tion would satisfy / = 1.0 and the MACHO mass mM = 1.0 Mq. 
The distributions arc not shown because they are essentially sim- 
ilar to Figure [5] except for the smaller integrated values (given in 
Tables EJand El . 



Table 6. Same as Table but this time using the MACHO 
collaboration's efficiency function B. 



Halo 


(Tcone ) tcxp 


(Tcone ) 


£cxp 




10~ 7 days 


10~ 6 


days 






events 








star -1 yr — 1 





#1 


9 1+2.2 
Z ' 1 -1.0 


84 


#2 


9 9+2.8 
z -1.0 


84 


#3 


l -°-0.7 


84 


#4 


1 9+ 21 

-0.9 


86 


#5 


1 9+ 13 


94 


#6 
#7 
#8 


1 S+ 1 - 9 
i-S-0.7 
i y+2.5 
L - -0.9 

i o+°- 8 


99 
89 
104 


mean 


1 7+ 1 ' 9 
i -'-0.8 


90 



0.7 
0.7 



+0.6 
0.3 
+0.8 
0.3 
+0.4 



0-5_o. 2 
6 +0 - 5 

u -°-0 .2 
6 +0 - 5 

U.D_ 2 

6 +0 - 6 
3 +0 - 2 



66 
71 
74 
74 
76 
71 
61 
81 



0.6 



+0.5 

0.2 



72 



Figure 1121 shows our computation for / and tum for 
every cone in all eight host haloes. The mean values are 
/ = 0.23^q' 13 and dim = 0.44io°i6> where the error limits 
show the mean scatter within a halo and contain 95 % of 
the values. 

Our values populate a pproximately the s ame areas as 
the contours in Figure 12 of lAlcock et alJ i2000T) with the ex- 
ception that our mass predictions do not extend to tum = 1.0 
Mq but stay mainly below tum = 0.6 Mq. Figure^] there- 
fore can be interpreted as (another) confirmation that down- 
scaled cluster sized CDM haloes acquired from cosmological 
iV-body simulations can be used in interpreting microlens- 
ing observations as well as analytical models — at least if 
the analysis is limited to the innermost regions. 

Note that the scatter in / is solely due to the variations 
in optical depth in different experiments, and thus, Figure 
1121 shows directly how optical depth varies between cones 
and even haloes. The scatter in myi contains additionally 
the variations in event rates between experiments. 

In the calculations above, we have assumed that r b s = 
3.38 x 10~ 8 is solely due to MACHOs in the dark halo of 
the Milky Way. However, if some of the observed events 
originate from elsewhere (e.g. the observations contain LMC 
self-lensing) , we have overestimated r b s . This leads to an 
overestimation of both / and tum, and in this sense, the 
values in Figure TTH are upper limits. 



6 DISCUSSION 

The benefits of using A-body haloes instead of analytical 
models are: (i) We do not have to make any initial assump- 
tions about the velocity distribution of the matter (other 
than limiting the circular velocity to 220 km/s), whereas 
analytical models have to adopt a Maxwellian distribution, 
(ii) We do not have to make any initial assumptions about 
the shape of the halo, whereas analytical models are always 
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-2.4 
-1.0 



90 
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educated guesses about the shape in form of some given pa- 
ram eters (e.g. triaxiality ; densit y profile, etc.). 

IWidrow fc Dubinskil l|l99Sft used a so called microlens- 
ing tube to get better number statistics for their event rates, 
and as they note, this "distorts the geometry of a realistic 
microlensing experiment" . However, our solution (the intro- 
duction of subparticles in accordance to the TSC mass as- 
sigment scheme) to the insufficient mass resolution of the 
cosmological simulation preserves the geometry. 

The downsides of using A-body haloes are: (i) We as- 
sume that the spatial and velocity distribution of MACHOs 
follows dark matter, (ii) The clusters we use are not as old 
as the Milky Way. 

Our models are based upon pure dark matter simula- 
tions which may not be appropriate if a significant fraction 
of the dark matter is composed of MACHOs, since they are 
composed of baryons, and baryonic physics has been ex- 
plicitly ignored! We thus implicitly assume that the spa- 
tial and velocity distribution of MACHOs follows the re- 
spective distributions of the underlying dark matter. Based 
upon these assumptions dark satellites comprised of MA- 
CHOs would be detectable through excess optical depth 
values and event duration anomalies in microlensing exper- 
iments. Simulations suggest that the majority of these dark 
satellites can be located as far as 400 kpc from the halo cen- 
tre. Thus, multiple background sources at distances over 400 
kpc would be needed to detect possible dark subhaloes in a 
Milky Way sized dark halo. For example, POINT -AGAPE 
jBelokurov et aUl2005L ICalchi Novati et alJl200rJ) is a sur- 
vey that is in principle able to detect even a dark MACHO 
satellite, on top of "free" MACHOs in the dark halo, because 
its source, M31, is distant enough. 

From the results, we can see that Halo #8 is too young 
to be used as a Milky Way dark halo model, but all the 
other haloes seem to behave well. The fact that Halo #8 
consists of three merging smaller haloes explains the peculiar 
values it produces. The other haloes are not experiencing 
any violent dynamical changes — a requirement we would 
expect a model of the Milky Way dark halo to fulfill. 

Obviously, our MACHO mass function is a <5-function. 
A more complex function would force us to calculate event 
durations and rates for individual MACHOs instead of sub- 
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particles. In this sense, subparticle values are only averages 
and some of the variation in the observables is lost. 

We chose not to use a complex MACHO mass function 
in this study because we wanted to concentrate on the halo 
structure effects. The variations of Y and t would be larger 
if we would use a MACHO mass function covering a large 
range of mass values. Thus, our results are as conservative 
as possible, and in real experiments even larger variations of 
event rates and durations could be expected. 

We recognize the fact that in the MACHO-lens scenario 
some of the lensing events towards the LMC could well be 
assigned to MACHOs within the dark halo of the LMC it- 
self. We have excluded the LMC's MACHO population in 
this study, but intend to investigate the characteristics and 
implications of such lens population in a separate, future 
study. Our plan is to model the dark matter halo of the 
LMC with some of the M ~ 10 9 Mq subhalos found near 
the cores of our iV-body haloes. 



7 CONCLUSIONS 

The purpose of this paper is to investigate the main mi- 
crolensing characteristics of JV-body dark matter haloes, ex- 
tracted from cosmological simulations and downscaled in 
size and mass to represent the dark halo of the Milky Way. 

We argue that analytical halo models are too simpli- 
fied in a number of ways in the case of microlensing where 
internal structures (in density or velocity distributions) or 
irregular halo shapes can alter the observed values signifi- 
cantly. 

We find that in general observables behave as expected 
from analytical models, resulting in fairly consistent values 
with seven out of the eight haloes with the "outsider" be- 
ing exceptionally young (3.42 Gyr) and undergoing a major 
merger of three smaller entties. As a result, this particular 
halo shows irregular behaviour in all tests and can not be 
considered as a valid model for the Milky Way dark halo. 

When individual experiments are examined in detail, 
we find that triaxiality and substructures can have a large 
effect on r and T. In some haloes, triaxiality can change the 
observed values by a factor as large as three. Substructure 
within haloes can also change r and T by a factor of two and 
furthermore reshape the event duration distribution notably. 

We also use our simulated values together with the MA- 
CHO collaboration's observations to find a preferred halo 
MACHO fraction (/) and individual MACHO mass (m M ). 
Our results are similar to the MACHO collaboration's own 
analysis where they used several different analytical models. 
In our analysis of / and thm, the scatter between different 
microlensing experiments is mainly due to triaxiality and no 
clear signs of frequent substructure signals can be seen. 
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Figure 9. Differential r and T values as a function of event duration. The largest Eigcn-value axis of the halo is along the z-axis, and 
the observers (40) are distributed uniformly onto the "solar sphere". Because each observer has 40 sources, the mean and the standard 
deviation are calculated from a total of 40 X 40 = 1600 cones. This is done by binning the subparticlc r and Y values for each cone into 
200 event duration bins and then calculating the mean and the standard deviation for each bin. Standard deviations are large due to 
triaxiality which causes the main differences in the mass and mass and velocity distributions between sightlines (cones). However, notice 
the Halo #1 standard deviation peak. It is due to a subhalo with a mass of 5.27 X 10 9 M (see text and Table^J. 
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Figure 10. The clear effect of triaxiality shown by r con c fractile boundaries which include 65 % of the values. The T CO nc values are 
binned in 25 bins by a, which is the angle between the source and one of the coordinate axes. Fractiles are calculated from t coi1c values 
within each bin. The major Eigen-value axis of the halo is aligned with the z-axis. Notice how T conc values grow when the source gets 
closer to the z-axis. 
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Figure 11. The same as Figure ITol but with T instead of t. 
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15 




Figure 12. The predicted MACHO mass (mivi) and MACHO halo fraction (f) values calculated for each cone. The predicted values 
are calculated as follows: / = T b s /r cono and m-M = (/r c0 ne/r o b s ) 2 . We derive the observational values from the observed event 
durations (criteria A, corrected only for blending, known populations JV e xp(-<4) = 2.67 subtracted). This leads to T b s = 3.38 X 10~ 8 
and r o b s = 1.69 X 10~ 7 events star - 1 yr -1 . The efficiency function is applied to the simulation values instead of the observed event 
durations. Our predicted values of / and mM are similar to the MACHO collaboration predictions. However, note the differences in 
the scatter area shapes between haloes. These differences are due to the fact that the haloes have different triaxial shapes. Moreover, 
substructure can cause a small number of individual data points to cluster together. 
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